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PURPOSE:  This  Coastal  and  Hydraulics  Engineering  Technical  Note  (CHETN)  provides  a 
literature  review  of  wave  dissipation  by  vegetation. 

INTRODUCTION:  Flooding  resulting  from  hurricanes  and  other  extreme  stonn  events  is  a 
prominent  risk  along  the  coasts.  These  coastal  areas  are  typically  of  low  elevation  and  relief, 
making  land  and  infrastructure  highly  susceptible  to  inundation  by  storm  surge  and  waves.  The 
severity  of  this  threat  is  exacerbated  by  sea  level  rise  and  a  possible  increase  in  stonn  frequency 
and  strength  due  to  climate  change.  Although  hard  protection  structures  such  as  levees  and 
floodwalls  reduce  flood  risk,  these  structures  may  fail  when  storm  conditions  exceed  the  design 
threshold.  There  is  a  general  consensus  that  wetlands,  which  often  serve  as  transition  zones 
between  open  water  and  dry  land,  could  act  as  buffers  and  reduce  storm  surge  and  propagating 
waves  substantially  before  they  encounter  coastal  development.  Unfortunately,  the  capability  of 
wetlands  to  serve  as  protection  during  extreme  storms  is  not  understood  fully  or  well 
documented;  furthermore,  water  level  and  wave  height  reductions  by  vegetation  are  studied  only 
in  low-energy  environments.  Nonetheless,  these  studies  present  methods  to  quantify  vegetation- 
induced  wave  attenuation  for  both  modeling  and  design.  This  technical  note  focuses  on  the 
damping  of  propagating  water  waves  by  vegetation,  but  also  discusses  surge  reduction  briefly. 
Although  waves  may  be  encountered  in  freshwater  environments  (e.g.,  boat  wakes,  lake  fetch, 
flood  waves,  etc.),  this  review  focuses  on  coastal  vegetation  and  resultant  effects  on  flood  and 
storm  damage  reduction. 

Water  waves  propagating  through  submerged  and  emergent  vegetation  lose  energy  by 
performing  work  on  the  vegetation  stems,  which  directly  results  in  smaller  wave  heights 
(Dalrymple  et  al.  1984).  Wave  attenuation  by  vegetation  is  a  function  of  vegetation  characteristics 
such  as  geometry,  buoyancy,  density,  stiffness,  and  spatial  coverage  as  well  as  wave  conditions 
such  as  incident  wave  height,  period,  and  direction.  Vegetation- wave  interactions  are  highly 
dynamic  in  that  the  vegetation  field  is  exposed  to  variable  wave  forcing  that  changes  with  time  as 
stems  bend,  flatten  to  the  bed,  or  are  washed  out.  As  evidenced  by  these  many  dependencies  and 
the  extensive  variety  of  coastal  plants,  the  variability  of  wave  damping  by  vegetation  is  large 
(Mendez  and  Losada  2004). 

Dean  and  Bender  (2006)  applied  linear  wave  theory  in  the  shallow-water  limit  to  show  that  wave 
setup,  a  component  of  storm  surge  resulting  from  the  transfer  of  breaking  wave  momentum  to  the 
water  column,  is  reduced  by  two-thirds  in  the  presence  of  vegetation  relative  to  conditions 
without  vegetation.  During  Hurricanes  Charley  and  Wilma,  water  levels  recorded  in  two  Florida 
mangrove  ecosystems  were  reduced  by  as  much  as  9.4  cm  per  km  inland.  Although  water  levels 
were  reduced  as  the  surge  moved  through  the  coastal  mangroves,  the  relative  contribution  of 
mangroves  was  still  unclear  (Krauss  et  al.  2009).  Numerical  simulations  by  Loder  et  al.  (2009)  and 
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Wamsley  et  al.  (2009)  lend  confidence  to  the  potential  of  vegetation  to  reduce  inundation  on  a 
landscape  scale.  However,  model  advancements  are  necessary  since  the  influences  of  vegetation 
are  often  approximated  as  bottom  friction. 

EXPERIMENTAL  RESULTS  OF  WAVE  DISSIPATION  BY  VEGETATION: 

Field  Studies.  The  ability  of  coastal  vegetation  to  dissipate  wave  energy  is  documented  and 
quantified  in  numerous  field  studies  focusing  on  low-energy  environments.  Field  studies  include 
those  conducted  along  intertidal  zones  of  England  (Moller  2006,  Cooper  2005,  Moller  and  Spencer 
2002,  Moller  et  al.1999),  coastal  mangrove  forests  of  Vietnam  (Quartel  et  al.  2007,  Mazda  et  al.  2006, 
Mazda  et  al.  1997),  salt  marshes  of  the  United  States  (Bradley  and  Houser  2009,  Knutson  et  al.  1982, 
Wayne  1976),  a  shallow  lake  in  Sweden  (Lovstedt  and  Larson  2010),  and  many  others.  Table  1 
summarizes  the  dominant  coastal  plant  species  and  the  wave  dissipation  observed  in  these  field 
experiments.  High  variability  in  wave  attenuation  is  recorded  in  these  limited  studies  as  a  result 
of  different  plant  species,  coverage,  and  wave  conditions. 


Table  1.  Field  studies  of  wave  attenuation  over  vegetation. 

Reference 

Transect  Length 
(m) 

Dominant  Plant  Species 

Average  Wave  Reduction 
(%  per  m) 

Wayne  (1976) 

20 

Spartina  alterniflora 

3.6 

20 

Thalassia  testudinum 

2.1 

Knutson  et  al.  (1982) 

30 

Spartina  alterniflora 

3.1 

Moller  et  al.  (1999) 

180 

Limonium  vulgare,  Aster  Tripolium,  Atriplex 
portulacoides,  Salicornia  spp.,  Spartina  spp., 
Suaeda  maritime,  Plantago  maritime, 
Puccinellia  maritima 

0.34 

Moller  and  Spencer  (2002) 

163 

Aster,  Suaeda,  Puccinellia,  Salicornia, 
Limonium  spp. 

0.54 

10 

Aster,  Suaeda,  Puccinellia,  Salicornia, 
Limonium  spp. 

4.38 

Cooper  (2005) 

300 

Puccinellia  maritima,  Salicorinia  europaea 

0.30 

250 

Atriplex  portulacoides,  Spartina  alterniflora 

0.26 

110 

Atriplex  portulacoides,  Salicorinia  europaea 

0.71 

Moller  (2006) 

10 

Spartina  anglica,  Salicornia  spp. 

1.8 

10 

Spartina  anglica,  Salicornia  spp. 

1.4 

10 

Salicornia  spp. 

1.0 

Quartel  et  al.  (2007) 

100 

Kandelia  candel,  Sonneratia  sp.,  Avicennia 
marina 

0.74 

Bradley  and  Houser  (2009) 

39 

Thalassia  testudinum 

0.77 

Lovstedt  and  Larson  (2010) 

over  first  5-14  m 
of  vegetation 

Phragmites  australis 

4.0 -5.0 

While  earlier  studies  tend  to  quantify  wave  dissipation  only  within  the  vegetation  (Knutson  et  al. 
1982,  Wayne  1976),  later  studies  compare  vegetated  and  unvegetated  areas  to  emphasize  the 
importance  of  maintaining  and  creating  wetlands  for  coastal  protection  (Quartel  et  al.  2007,  Mazda 
et  al.  2006,  Cooper  2005,  Moller  and  Spencer  2002,  Moller  et  al.  1999).  A  year-long  study  by  Cooper 
(2005)  found  that  wave  attenuation  over  salt  marshes  was  at  least  two  times  greater  than  over 
mudflats.  The  largest  discrepancy  along  one  transect  was  a  wave  height  reduction  of  0.30%  per 
m  across  a  salt  marsh  versus  0.03%  per  m  across  a  mudflat.  Attenuation  over  another  UK  salt 
marsh  was  approximately  50%  greater  relative  to  a  mudflat  considering  the  same  water  depths 
(Moller  et  al.  1999).  Coastal  mangroves  in  the  Red  River  Delta,  Vietnam,  were  found  to  reduce 
wave  heights  5-7.5  times  more  than  bottom  friction  alone  (Quartel  et  al.  2007).  The  other 
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mentioned  citations  reflect  similar  results  where  attenuation  was  considerably  greater  over 
vegetation  than  bare  transects. 

Details  of  the  pattern  of  wave  attenuation  as  well  as  the  behavior  of  the  vegetation  are  of  interest. 
Wave  heights  were  observed  to  decay  exponentially  through  vegetation  beds  with  the  most  rapid 
reduction  occurring  within  the  first  few  meters  of  permanent  vegetation  cover.  Knutson  et  al. 
(1982)  reported  a  20%  wave  height  reduction  per  m  within  the  first  2.5  m  while  Moller  and 
Spencer  (2002)  reported  values  of  1 . 14%  and  2.12%  per  m  within  the  first  10  m  for  two  transects. 

A  recent  study  by  Bradley  and  Houser  (2009)  characterized  the  relative  movement  of  seagrass  blades 
in  oscillatory  flow  to  address  the  importance  of  blade  movement  to  wave  attenuation.  Wave 
attenuation  tended  to  decrease  as  incident  wave  height  increased  and  conditions  became  more 
turbulent  (200  <  Rea  <  800  where  Rea  is  the  stem  Reynolds  number  defined  as  Red=  ucd! v  and  uc  is 
the  characteristic  velocity  acting  on  the  vegetation  -  typically  the  maximum  horizontal  velocity,  d 
is  the  plant  width  nonnal  to  the  flow,  and  v  is  kinematic  viscosity).  Video  observations  showed 
seagrass  blades  swaying  over  the  entire  wave  cycle  when  exposed  to  lower  orbital  velocities. 
However,  at  higher  orbital  velocities,  the  blades  extended  in  the  direction  of  flow  for  the  majority 
of  the  cycle  and  became  streamlined,  leading  to  reduced  drag  and  less  attenuation.  The  blades  were 
also  found  to  move  in-phase  with  waves  at  a  lower  frequency  (0.38  Hz)  but  out-of-phase  with 
waves  at  the  higher  peak  frequency  (0.67  Hz).  This  behavior  suggests  higher  frequencies  (shorter 
wave  periods)  in  a  spectrum  should  experience  more  attenuation  by  seagrass;  these  results  are 
supported  by  Lowe  et  al.  (2007)  who  found  that  shorter-period  waves  penetrated  more  readily  into  a 
canopy  of  staggered  plastic  cylinders  and  were  dissipated  more  effectively. 

Lovstedt  and  Larson  (2010)  studied  wave  dissipation  by  reeds  in  a  shallow  lake.  They  analyzed 
changes  in  the  wave  height  probability  density  function  (PDF)  for  waves  traveling  through 
vegetation.  It  was  expected  that  the  Rayleigh  distribution,  which  is  often  used  to  characterize 
random  variation  in  wave  height,  would  be  transfonned  as  higher  waves  were  preferentially 
damped.  Although  a  slight  shift  was  observed  toward  the  transfonned  PDF,  especially  in  shallow 
water  and  for  longer  distances  into  the  vegetation,  the  shift  was  not  significant  and  merits 
additional  attention. 

Laboratory  Studies.  Vegetation- wave  interactions  have  also  been  studied  under  controlled 
conditions  in  laboratory  flume  studies  with  either  natural  vegetation  (Tschirky  et  al.  2000,  Fonseca 
and  Cahalan  1992)  or  artificial  vegetation  (Cavallaro  et  al.  2010,  Augustin  et  al.  2009,  Lima  et  al.  2006, 
Lovas  and  Torum  2000,  Dubi  and  Torum  1996)  simulated  by  various  elements.  The  “vegetation” 
selected  for  each  experiment  is  shown  in  Table  2.  These  flume  experiments  often  isolate  the 
influence  of  a  single  plant  or  wave  property  on  wave  dissipation  by  varying  one  parameter  of 
interest  while  holding  other  parameters  constant.  Tschirky  (2000)  identified  several  trends  with 
respect  to  wave  attenuation  through  wetland  plants: 

1 .  total  wave  attenuation  increased  with  longer  plant  beds, 

2.  higher  stem  densities  resulted  in  greater  attenuation, 

3.  attenuation  slightly  increased  with  larger  incident  wave  heights, 

4.  attenuation  varied  inversely  with  water  depth,  and 

5.  no  clear  trend  with  respect  to  wave  period  was  observed. 
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Table  2.  Stem  elements  for  flume  experiments. 

Reference 

Element 

Cavallaro  et  al.  (2010) 

polyethylene  Posidonia  oceanica 

Augustin  et  al.  (2009) 

cylindrical  wooden  dowels  for  Spartina  altemiflora 

polyethylene  foam  tubing  for  Spartina  altemiflora 

Lima  et  al.  (2006) 

flexible  nylon  rope  for  Brachiaria  subquadripara 

Lovas  and  Torum  (2000) 

plastic  Laminaria  hyperborea 

Tschirky  et  al.  (2000) 

harvested  Scirupus  americanus 

Dubi  and  Torum  (1996) 

plastic  Laminaria  hyperborea 

Fonseca  and  Cahalan  (1992) 

harvested  Halodule  wrightii,  Syringodium  filiforme,  Thalassia  testudinum, 
Zostera  marina 

Parametric  analyses  of  wave  attenuation  through  vegetation  have  been  completed  by  a  number  of 
other  authors  with  several  supporting  Tschirky’s  (2000)  conclusions.  Lima  et  al.  (2006)  reported 
attenuation  values  of  0.83-1.67%  per  m  for  400  stems/nr  versus  5%  per  m  for  600  stems/m". 
Unlike  the  effect  of  stem  density,  the  influence  of  wave  height  is  less  clear.  The  preliminary 
results  of  Cavallaro  et  al.  (2010)  also  exhibit  an  increase  in  wave  attenuation  with  larger  incident 
wave  height;  however,  Bradley  and  Houser  (2009)  observed  less  wave  attenuation  with  increasing 
wave  height.  Lovas  and  Torum  (2000)  found  that  although  the  smallest  maximum  wave  height 
values  were  reduced  the  most,  the  largest  difference  between  unvegetated  and  vegetated 
conditions  was  observed  for  the  largest  offshore  wave  heights. 

A  convenient  way  to  address  the  single  effect  of  water  depth  on  wave  attenuation  is  to  define  the 
relative  stem  length.  The  relative  stem  length  represents  the  percentage  of  the  water  column 
occupied  by  the  stem  and  is  the  ratio  of  stem  length  (/,)  to  water  depth  (h).  Augustin  et  al.  (2009), 
Dubi  and  Torum  (1996),  and  Fonseca  and  Cahalan  (1992)  observed  substantially  more  wave 
reduction  under  emergent  conditions  ( IJh  >  1.0),  with  dissipation  decreasing  as  water  depth 
exceeded  canopy  height  {IJh  <  1.0).  Based  on  theoretical  and  experimental  results,  Dubi  and 
Torum  (1996)  propose  a  50%  reduction  in  significant  wave  heights  over  74,  216,  495,  and  1390  m 
for  IJh  =  0.50,  0.33,  0.25,  and  0.2,  respectively,  for  a  stem  density  of  12  stems/m".  Augustin  et  al. 
(2009)  observed  50-200%  greater  wave  attenuation  per  wave  length  over  emergent  than  near- 
emergent  ( IJh  =  0.75)  stems.  Field  results  of  wave  attenuation  through  mangroves  reveal  a 
slightly  different  pattern  than  that  of  seagrasses.  Dissipation  decreased  as  water  exceeded  the 
height  of  the  pneumatophores  (an  extensive,  hardy  root  system),  similar  to  when  water  depth 
exceeded  the  height  of  the  seagrass  canopy.  However,  when  the  water  level  reached  the  height  of 
the  mangrove  branches  and  leaves,  the  rate  of  wave  reduction  increased  again  (Mazda  et  al.  2006). 
Clearly  the  amount  of  wave  dissipation  is  dependent  on  the  plant  structure  (relative  to  the  water 
depth).  The  value  of  4  may  change  over  time  as  the  water  depth  varies  (due  to  tide  or  stonn 
surge)  or  the  vegetation  flexes. 

A  relationship  between  relative  vegetation  height  and  wave  dissipation  is  expected  given  that 
wave  horizontal  particle  velocities,  u(z),  are  highest  near  the  crest  of  the  wave  with  velocities 
decreasing  towards  the  bed  (depending  on  water  depth  and  wave  period)  (Dean  and  Dalrymple 
1991).  The  highest  velocities  are  impeded  as  stems  approach  the  surface  of  the  water,  leading  to 
greater  drag  and  wave  dissipation. 

Intuitively,  the  impact  of  vegetation  is  greatest  during  emergent  conditions  where  the  entire 
velocity  profile  is  impeded.  As  stems  become  submerged,  the  highest  orbital  velocities  pass  over 
the  canopy.  Only  slower  wave  velocities  occurring  over  the  lower  part  of  the  water  column  are 
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impeded,  thus,  less  drag  is  generated,  which  translates  to  less  wave  attenuation.  An  illustration  of 
emergent  and  submerged  conditions  is  presented  in  Figure  1.  Wetlands  may  be  less  effective  in 
attenuating  waves  during  stonn  events  due  to  submergence,  but  Fonseca  and  Cahalan  (1992) 
suggest  that  even  a  relatively  low  rate  of  wave  reduction  across  wetlands  can  be  substantial  when 
applied  over  long  distances. 


Figure  1.  Interactions  of  particle  velocities  with  emergent  (a)  and  submerged  (b)  stems. 


In  comparison  to  the  other  parameters,  the  effect  of  wave  period  on  dissipation  is  the  most 
contradictory  and  unresolved.  Moller  et  al.  (1999)  observed  salt  marshes  reducing  wave  energy 
at  all  frequencies  by  the  same  degree.  However,  other  field  and  laboratory  studies  indicate  a 
dominating  frequency-dependent  component  to  wave  attenuation  over  vegetation.  As  previously 
mentioned,  Bradley  and  Houser  (2009)  and  Lowe  et  al.  (2007)  both  describe  seagrass  as  a  low- 
pass  filter  in  which  higher  frequency  components  are  dissipated  more  effectively.  Cavallaro  et  al. 
(2010)  and  Fonseca  and  Cahalan  (1992)  address  wave  period  by  defining  the  relative  water  depth 
(h/L)  where  L  is  wave  length.  This  representation  allows  the  investigation  of  wave  period  as  a 
function  of  water  depth,  which  determines  whether  a  wave  is  shallow  water  (h/L  <  0.05), 
intermediate  (0.05  <  h/L  <  0.5),  or  deepwater  (h/L  >  0.5)  (Sorensen  2006).  Wave  attenuation 
decreased  as  waves  approached  the  deepwater  regime  for  both  Cavallaro  et  al.  (2010)  and 
Fonseca  and  Cahalan  (1992). 

Shallow-water  and  deepwater  waves  have  different  particle  trajectories,  and  these  trajectories 
affect  the  amount  of  drag  generated  and  wave  energy  dissipated.  Diagrams  of  particle  trajectories 
for  shallow-water  and  deepwater  waves  are  shown  in  Figure  2.  The  particle  path  for  a  deepwater 
wave  is  a  circular  orbit  that  decays  exponentially  with  negligible  particle  motion  at  an  elevation 
below  the  surface  >  -0.5Z.  The  horizontal  excursion  for  a  shallow-water  wave  is  constant  over 
depth  but  the  vertical  excursion  decreases  linearly  from  the  surface  to  zero  at  the  bottom,  resulting 
in  an  elliptic  orbit  (Sorenson  2006,  Dean  and  Dalrymple  1991).  The  interaction  of  waves  with 
vegetation  can  vary  with  relative  depth  (or  wave  period),  and  it  is  anticipated  that  vegetation  will 
affect  shallow-water  waves  more  than  deepwater  waves  because  shallow  waves  tend  to  have  larger 
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velocities  throughout  the  water  column  and  near  the  bed.  It  is  important  to  note  that  waves  interact 
with  the  bottom  based  on  relative  depth,  so  waves  with  longer  wave  lengths  become  shallow-water 
waves  and  interact  with  the  bottom  at  deeper  depths  (Koch  et  al.  2006). 


MODELING  WAVE-VEGETATION  INTERACTIONS:  Various  models  and  model  extensions 
for  wave-vegetation  interactions  have  been  proposed.  While  some  models  approximate  vegetation 
with  higher  bottom  friction  factors  (Camlield  1977,  Moller  et  al.  1999),  the  majority  implement  an 
empirical  drag  coefficient  to  estimate  the  wave-induced  drag  forces  along  the  plant  stem 
(Dalrymple  et  al.1984,  Asano  et  al.1992,  Kobayashi  et  al.  1 993,  Mendez  and  Losada  1999,  Mendez 
and  Losada  2004,  Lima  et  al.  2006,  Lowe  et  al.  2007,  Mullamey  et  al.  2010).  Massel  (1999)  and 
Vo-Luong  (2008)  focus  exclusively  on  modeling  the  interactions  between  propagating  water 
waves  and  mangroves.  Applying  these  models  to  appropriate  physical  conditions,  based  on  the 
biomechanics  of  the  plant  species,  is  crucial.  For  example,  reed  plants  such  as  S.  alterniflora  can  be 
simulated  as  rigid,  vertical  cylinders  whereas  flexible  vegetation  such  as  kelp  requires  a  more 
complex  drag  fonnulation.  Figure  3  is  a  diagram  of  typically  used  parameters  in  wave-vegetation 
modeling,  where  As  is  average  stem  spacing,  d  is  average  stem  diameter,  4  is  average  stem  length, 
N  is  the  average  number  of  stems  per  unit  area  (average  stem  density,  stems/m2),  Fx  is  total 
horizontal  force  per  unit  volume  on  a  stem  array,  and  h  is  water  depth. 


Figure  3.  Vegetation  modeling  parameters  (modified  from  Dean 
and  Bender  2006). 
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One  of  the  first  hydrodynamic  models  for  wave  attenuation  was  proposed  by  Dean  (1979)  based 
on  empirical  estimates  of  fluid  drag  forces  acting  on  vertical,  rigid  cylinders.  The  model  for  the 
damping  of  incident  wave  height  ( H0 )  by  coastal  plants  was  the  following: 

H  _  1 

JT0~  1  +  A'x 


where 


A' 


CDd 
6kA  s2h 


H 


o 


H  is  the  local  wave  height,  x  is  the  cross-shore  distance  waves  propagate  through  the  vegetation 
field  (Figure  3),  and  A  ’  is  the  rate  of  wave  energy  dissipation.  This  form  of  equation  (H/I I a)  is 
used  by  several  authors  with  differing  dissipation  rate  formulations.  The  bulk  or  average  drag 
coefficient  Cd  for  the  entire  plant  field  was  assumed  to  be  constant  over  the  depth  and  was 
approximated  as  1 .0  to  describe  the  forces  associated  with  smooth,  rigid  cylinders.  Knutson  et  al. 
(1982)  recognized  that  the  Dean  model  did  not  describe  the  responses  of  plants  to  wave  forcing, 
such  as  swaying,  and  incorporated  an  empirical  plant  drag  coefficient,  Cp,  to  account  for  the 
difference  between  rigid  cylinders  and  plants, 

H  _  1 

JT0~1  +  A^ 

where 


A  =  — e—^—Ho 
QkAs21i 

A  calibrated  value  of  Cp  =  5  minimized  the  error  between  predicted  and  measured  wave  heights 
for  data  collected  in  S.  alterniflora  marshes  of  Chesapeake  Bay. 

Dalrymple  et  al.  (1984)  formulated  an  algebraic  dissipation  equation  using  linear  wave  theory  and 
conservation  of  wave  energy  by  approximating  a  vegetation  bed  as  an  array  of  rigid,  vertical 
cylinders.  Wave  dissipation  is  assumed  to  only  be  a  function  of  the  horizontal  drag  force  per  unit 
volume  (inertial  component  neglected),  which  is  summed  over  all  the  stems  and  is  expressed  as, 

1  .  . 

Fx  —  —pCDNdu  \u\ 

where  u  is  the  horizontal  fluid  velocity  due  to  wave  motion  and  p  is  the  fluid  density.  The  time- 
averaged  energy  dissipation  resulting  from  vegetation  is  the  product  of  Fx  and  u,  and  is 
integrated  over  the  stem  length.  The  solution  for  wave  height  in  a  vegetation  field  was  given  by: 

H  _  1 

Ho  1  -\-ax 
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where 


a  =  — CD  Ad  (sinh  3  kls  +  3  sinh  kls ) 
3  TZ 


4k 


3  sinh  A7?  (sinh  2  kh  +  2kh) 


H, 


and  k  is  the  wave  number.  The  model  accounts  for  arbitrary  water  depth  and  stem  length, 
allowing  for  both  submerged  and  emergent  vegetation,  and  the  bulk  drag  coefficient  CD  is 
empirically  detennined  to  account  for  varying  responses  of  different  plant  species  to  wave 
forcing  (i.e.,  calibrated  to  different  plant  structure  or  flexibility).  Mendez  and  Losada  (2004) 
expanded  upon  Dalrymple  et  al.  (1984)  and  developed  an  empirical  model  to  estimate 
monochromatic  and  random  wave  transfonnations  over  mildly-sloping  vegetation  fields  under 
nonbreaking  and  breaking  conditions.  The  random  wave  expansions  were  validated  for  the 
artificial  Laminaria  hyperborea  kelp  experiments  of  Dubi  (1995)  and  Lovas  (2000)  for 
nonbreaking  and  breaking  waves,  respectively. 

Kobayashi  et  al.  (1993)  derived  an  analytical  solution  based  on  the  continuity  and  linearized 
momentum  equations  for  small-amplitude  monochromatic  wave  attenuation  over  submerged 
vegetation.  The  continuity  equation  within  the  vegetation  was  given  by, 


du 

dx 


dw 

Ik 


0 


and  the  linearized  momentum  equations  per  unit  water  volume  was  expressed  as, 


du  _  dp 
p  dt  dx 


dw  dp  ^ 

p - = - F_ 

dt  dz  z 

where  w  is  vertical  particle  velocity,  z  is  the  vertical  coordinate,  t  is  time,  Fz  is  the  total  vertical 
force  per  unit  volume,  and  p  is  dynamic  pressure.  Like  Dalrymple  et  al.  (1984),  plant  stems  are 
approximated  as  vertical  cylinders,  and  Fz  is  assumed  negligible  (Fz  ~  0)  when  compared  to  Fx. 
However,  while  Dalrymple  et  al.  (1984)  implements  a  quadratic  drag  formulation,  Kobayashi  et 
al.  (1993)  linearizes  the  horizontal  drag  force  as  a  function  of  fluid  particle  velocity.  The  local 
wave  height  is  assumed  to  decay  exponentially  with  propagation  through  a  vegetation  bed 
according  to  the  following  fonn: 


=  exp(— k.x) 


where  A,  is  the  exponential  damping  coefficient.  Calibrated  values  of  Cd  are  obtained  by  matching 
calculated  values  of  the  nonnalized  decay  coefficient,  k{  =  ki  h,  to  measured  k{  values.  The 
exponential  decay  solution  captured  adequately  the  wave  trends  observed  by  Asano  et  al.  (1988). 


8 


ERDC/CHL  CHETN-l-82 
September  201 1 


Asano  et  al.  (1992)  refined  the  analytical  solution  developed  by  Kobayashi  et  al.  (1993)  to 
include  small  amplitude  vegetation  motion.  A  more  correct  estimation  of  Fx  is  calculated  using 
the  relative  velocity  (ur)  between  the  horizontal  particle  velocity  and  swaying  velocity  of  the 
stem  ( uv ), 


Fx  =  ~pCDNdur  | ur 


where  ur  =  (u  -  uv).  The  drag  force  is  then  linearized  as  in  Kobayashi  et  al.  (1993).  Each 
individual  stem  was  modeled  as  a  cantilever  beam,  fixed  at  the  bottom,  with  swaying  modeled  as 
a  forced  vibration  with  one  degree  of  freedom  where  buoyancy  and  stiffness  are  the  restoring 
forces.  The  solution  of  this  problem  is  determined  iteratively  and  indicated  by  a  converged 
solution.  Comparisons  between  the  present  model  and  the  results  of  Asano  et  al.  (1988)  yielded  a 
better  agreement  than  the  solution  of  Kobayashi  et  al.  (1993).  The  wave  decay  solutions  of 
Kobayashi  et  al.  (1993)  and  Asano  et  al.  (1992)  are  further  extended  by  Mendez  and  Losada 
(1999).  The  complete  wave  field  is  included  by  separating  the  system  into  four  regions  (region 
seaward  of  the  vegetation,  region  above  the  vegetation,  region  behind  the  vegetation,  and  within 
the  vegetation)  and  obtaining  a  velocity  potential  for  each  region  by  solving  the  boundary  value 
problem.  This  method  accounts  for  reflection  and  transmission  of  waves  as  they  encounter  the 
interfaces  between  the  seaward,  vegetation,  and  leeward  regions.  Mendez  and  Losada  (1999) 
describe  Fx  in  its  entirety  using  a  linear  summation  of  the  drag  force  and  inertial  force,  which  is 
caused  by  fluid  and  vegetation  accelerations, 


Fx  =  -pCDNdur  | ur  |  +  pNdtv 


Cdu 


dt 


\cm 


i) 


d2c 

dt2 


where  tv  is  the  stem  width,  Cm  is  the  inertia  coefficient,  and  C  is  the  horizontal  displacement  of 
the  vegetation  stand  with  respect  to  the  vertical  axis.  The  model  solution  compared  well  with  the 
experimental  results  of  Asano  et  al.  (1988)  and  captured  accurately  the  modulation  in  wave 
heights  resulting  from  reflection  at  the  vegetation  interfaces.  The  formulation  was  then  extended 
to  irregular  waves,  and  model  results  compared  well  to  Dubi  (1995). 

Lima  et  al.  (2006)  formulated  an  algebraic  model  for  monochromatic  wave  damping  over  highly 
flexible  vegetation  (nylon  string)  using  nonlinear  stream  function  wave  theory.  Unlike  the 
previously  discussed  models  where  Fx  is  simply  a  summation  of  the  drag  contributed  by 
individual  stems,  the  resistances  resulting  from  interactions  between  stems,  such  as  entangling,  is 
accounted  for  using  an  ensemble  interaction  coefficient  (r).  The  horizontal  drag  force  in  which 
all  stem  elements  are  identical  is  defined  as, 

Fx  =  tNF * 


where  Fx  is  the  force  exerted  by  a  single  stem  on  the  fluid.  Each  nylon  string  is  treated  as  a  string 

of  equally  spaced  nodes  with  freedom  of  movement  in  the  x-z  plane,  and  is  described  by  a  system 
of  three  equations;  a  constitutive  equation,  and  the  momentum  balances  in  the  horizontal  and 
vertical  directions.  The  solution  of  the  system  provides  the  connecting  forces  between  nodes  and 
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the  horizontal  and  vertical  displacement  of  each  node.  The  velocity  and  acceleration  of  individual 
stems  were  determined  from  video  and  allowed  Cd  to  be  best  fit  according  to  the  expression. 


FT  =  2pCDA*Ur  k 


where  Ft  is  the  hydrodynamic  forcing  and  As  is  the  area  of  the  stem  normal  to  the  relative 
velocity.  According  to  Newton’s  Law,  the  resisting  force  F'x  of  the  stem  against  the  flow  equals 

in  magnitude  the  force  Ft  exerted  by  the  flow  on  the  stem.  The  drag  force  for  a  single  element  is 
then  obtained  as  a  function  of  flow  and  plant  properties  by  dimensional  analysis, 


Prd 


F;  —  0.133- 


H a 


where  pv  is  the  density  of  the  vegetation  element  and  T  is  wave  period.  Following  calibration 
with  single  elements,  the  ensemble  coefficient  and  wave  decay  model  are  derived  from  the 
behavior  of  the  stem  field.  The  ensemble  coefficient  is  given  as, 


2  ,3/2 


r  =  7.44- 


L  d 


HJ2As 


1/2 


Substituting  r,  the  wave  height  decay  model  for  waves  traveling  through  a  vegetation  field  with 
flexible  elements  is  the  following: 


H  _  1 

IF Vx 


where 


P 


L2 


7.916 


5/2  A  75/4 


d°,zN- 


gT2(An  —  l)/2 


Pv 

P 


and  n  is  the  ratio  between  wave  group  speed  and  celerity  given  as  (Dean  and  Dalrymple  1991), 


sinh  2kh 

Lima  et  al.  (2006)  concluded  that  the  resistance  imposed  by  a  group  of  stems  was  on  average 
four  times  higher  than  summing  individual  forces  for  all  stems,  and  suggested  that  stem 
interactions  are  of  a  greater  importance  than  previously  thought. 

Mullamey  et  al.  (2010)  derived  an  analytical  solution  for  the  wave-forced  movement  of  single¬ 
stem  vegetation  and  compared  the  results  against  observed  motion  of  Schoenoplectus  americanus. 
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The  majority  of  measured  relationships  between  water  and  stem  motion  were  in  good  agreement 
with  model  theory.  Theoretical  analyses  identified  stiffness  as  an  important  parameter  influencing 
wave  attenuation.  Low  stiffness  stems  moved  with  the  surrounding  water;  in  contrast,  stiffer  stems 
led  the  surrounding  water  by  90°  and  have  little  mobility.  Mullamey  et  al.  (2010)  defined  stiffness 
using  vegetation  biomechanical  and  wave  properties;  low  stiffness  is  associated  with  a  low 
Young’s  modulus,  long  thin  stems,  and  high-energy,  high-frequency  waves.  Simulated  dissipation 
predicted  for  flexible  stems  was  approximately  30%  that  of  rigid  stems. 

VEGETATION  DRAG  COEFFICIENTS:  The  hydrodynamic  forces  acting  on  the  vegetation 
due  to  waves  is  F  =  (Fx,  0,  Fz),  where  it  is  usually  assumed  that  Fz  is  negligible  in  comparison  to  Fx 
(shallow  water).  The  horizontal  force  Fx  is  described  typically  using  a  Morison-based  equation. 
The  Morison  equation  describes  the  forces  induced  by  a  solid  body  in  oscillatory  flow  and  is 
comprised  of  two  parts:  a  drag  force  (FD)  and  an  inertial  force  (F/).  The  Morison  equation  is  given 
as, 


F  =  FD  +  F1=VpCDA  \u\u  +  pCmV -± 

where  Cd  is  the  drag  coefficient  for  a  singular  element,  A  is  the  cross-sectional  area  of  the  body 
perpendicular  to  the  flow,  and  V  is  the  body’s  volume  (Morison  et  al.  1950).  The  values  of  the 
drag  and  inertia  coefficients  are  determined  by  experiments.  In  the  case  of  a  cylinder  of  diameter 
d  is  oscillatory  flow,  the  reference  area  per  unit  length  is  A  =  d.  Only  the  drag  force  is  considered 
in  the  majority  of  the  above  models,  and  an  average  or  bulk  drag  coefficient  Cd  serves  as  the 
calibration  parameter  to  minimize  error  between  measured  and  predicted  wave  heights.  The  bulk 
drag  coefficient  allows  for  the  total  drag  force  on  a  vegetation  field  to  be  a  summation  of  the 
drag  forces  contributed  by  individual  stems. 

Defining  a  generalized  value  to  describe  all  plant-induced  drag  is  impossible  since  Cd  is  highly 
dependent  on  hydrodynamic  and  plant  characteristics.  However,  the  generalized  behavior  of  the 
vegetation  bulk  drag  coefficient  can  be  understood  by  examining  empirical  fonnulas.  These 
empirical  equations  attempt  to  predict  appropriate  Cd  values  for  specific  plant  types  by  defining 
Cd  as  a  function  of  nondimensional  flow  parameters. 

One  parameter  found  to  influence  the  bulk  drag  coefficient  is  the  stem  Reynolds  number, 


v 


To  calculate  the  stem  Reynolds  number,  Kobayashi  et  al.  (1993)  substituted  the  maximum 
horizontal  velocity  given  by  linear  theory  for  uc, 

_  gHk  cosh  k{h  +  z) 

Wmax  “  2 co  cosh  kh 

where  co  is  wave  angular  frequency  and  defined  as  co  =  (2jt/T).  Kobayashi  et  al.  (1993)  plotted 
CD  values  calibrated  for  the  artificial  kelp  experiments  of  Asano  et  al.  (1988)  against  the 
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corresponding  stem  Reynolds  number  and  found  the  following  relationship  for  2,200  <  Red  < 
18,000: 


CD  —  0.08  -|- 


/  \2.4 

2200 

,  Red  . 


The  bulk  drag  coefficient  decreases  with  increasing  Red  and  approaches  the  order  of  0. 1  for  large 
Red.  Mendez  et  al.  (1999),  who  extended  the  model  of  Kobayashi  et  al.  (1993)  to  include 
vegetation  motion  and  interfaces,  also  used  Asano  et  al.  (1988)  data  for  model  validation. 
However,  the  empirical  relationships  vary  from  that  proposed  by  Kobayashi  et  al.  (1993)  and 
differ  based  on  whether  vegetation  motion  is  considered, 

Swaying  neglected: 


C D  —  0.08  + 


\2.2 

2200 

,  Red  , 


for  200  <  Red  <  15,500 


Swaying  considered: 


CD  =  0.40  + 


\2.9 

4600 

,  Red  . 


for  2300  <Red  <20,000 


The  “no  swaying”  equation  proposed  by  Mendez  et  al.  (1999)  yielded  a  20%  improvement  in  the 
correlation  coefficient  over  that  of  Kobayashi  et  al.  (1993)  and  resulted  in  less  scatter.  These 
relationships  compared  with  Kobayashi  et  al.  (1993)  are  presented  in  Figure  4.  It  is  clear  that  the 
inclusion  of  plant  motion  results  in  a  larger  bulk  drag  coefficient  given  the  same  stem  Reynolds 
number.  This  occurs  because  the  characteristic  velocity  substituted  into  the  drag  fonnulation  is 
smaller  when  considering  plant  motion.  Kobayashi  et  al.  (1993)  neglected  swaying  and 
substituted  maximum  horizontal  particle  velocity  given  by  linear  wave  theory  while  Mendez  at 
al.  (1999)  substituted  the  relative  velocity  between  the  fluid  and  plant.  This  relative  velocity  is 
smaller  than  that  given  by  linear  wave  theory  and,  therefore,  a  higher  value  of  Cd  is  required  to 
obtain  the  same  amount  of  attenuation. 

Mendez  and  Losada  (2004)  attempted  to  parameterize  Cd  for  the  artificial  kelp  experiments  of 
Dubi  (1995)  as  a  function  of  the  local  Keulegan-Carpenter  number,  K,  defined  as, 


d 

The  maximum  velocity  uc  for  random  waves  was  defined  using  the  root-mean-square  wave 
height  Hrms  and  peak  wave  period  Tp.  However,  Mendez  and  Losada  (2004)  observed  scatter 
when  Cd  was  considered  to  be  a  function  of  K  alone  and  found  that  the  drag  coefficient  also 
depended  on  the  relative  stem  length.  To  allow  Cd  to  be  a  function  of  K  and  the  relative  stem 
length,  a  modified  Keulegan-Carpenter  number  (Q)  was  defined  and  the  relationship  between  CD 
and  Q  was  defined  as  the  following  for  7  <  Q  <  172: 
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Figure  4.  Empirical  relationships  between  Red  and  CD. 


_  exp(-O.O1380 

D  _  n  o 


where 


Q  =  —^~ 

(/,/*) 

The  bulk  drag  coefficient  decreases  with  increasing  A"  and  increases  with  larger  relative  stem 
lengths  given  the  same  Keulegan-Carpenter  number.  An  increase  in  the  drag  coefficient  with 
higher  relative  stem  length  is  expected  since  a  larger  percentage  of  the  water  column  is  occupied 
by  the  plant.  This  trend  can  be  seen  in  Figure  5.  Augustin  et  al.  (2009)  found  CD  was  better 
correlated  with  the  stem  Reynolds  number  during  emergent  conditions  and  the  Keulegan-Carpenter 
number  during  near-emergent  conditions.  However,  Bradley  and  Houser  (2009)  found  no 
improvement  in  parameterizing  Co  with  the  Keulegan-Carpenter  number  versus  the  stem  Reynolds 
number. 

Hydrodynamic  studies  with  steady  unidirectional  flows  and  vegetation  suggest  interactions 
between  cylinders  also  affect  the  bulk  drag  coefficient.  Nepf  (1999)  found  that  CD  decreased  as 
stem  density  increased  as  a  result  of  turbulence-scale  wake  sheltering.  Wake  sheltering  is  an 
interaction  between  upstream  and  downstream  cylinders  whereby  upstream  stem  wakes  reduce 
the  drag  on  downstream  stems.  Variations  in  vegetation  geometry  and  spatial  coverage  might 
also  change  the  attenuation  capacity  of  a  vegetation  field.  Future  controlled  studies  investigating 
how  stem  spatial  variation,  and  hence  turbulence  flow  fields  and  vortex  interactions,  relate  to 
wave  attenuation  are  recommended  Irish  et  al.  (2008). 
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Myrhaug  et  al.  (2009)  outlines  a  practical  method  for  estimating  the  drag  force  on  a  vegetation 
field  subjected  to  nonlinear  random  waves.  This  is  achieved  using  a  common  drag  formulation 
and  the  empirical  relationships  for  Cd  proposed  by  Mendez  et  al.  (1999).  Like  the  other 
formulations,  wave-induced  drag  forces  are  described  using  a  Morison-type  equation  with 
vegetation  motion  and  inertial  forces  neglected  where  the  maximum  horizontal  drag  force  within 
a  wave  cycle  is  given  by, 


^nax  =  ~pCDNdu 


2 

max 


The  relationships  proposed  by  Mendez  et  al.  (1999)  are  of  the  following  form: 

P 


CD  —  a  + 


Rc , 


where  for  rigid  plants, 


(a,  p,  y)  =  (0.08,  2200,  2.2)  for  200  <  Red  <  15,500 


and  for  swaying  plants, 

(a,  p,  y)  =  (0.40,  4600,  2.90)  for  2300  <  Red  <  20,000 
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Nonlinear  wave  effects  are  taken  into  account  using  Stokes  second-order  wave  theory.  Stochastic 
and  detenninistic  methods  to  determine  the  drag  coefficient  are  presented,  although  Myrhaug  et 
al.  (2009)  emphasize  using  the  stochastic  method.  An  estimate  for  the  root-mean-square  drag 
coefficient  is  the  following: 


C 


Drms 


=  rJ~B 


1  +  kpan 


sA 

rB, 


where  kp  is  the  wave  number  associated  with  the  spectral  peak  and  anm  is  the  root-mean-square 
wave  amplitude.  The  other  variables  (s,  r,  A,  B,  Rrms)  are  given  as, 


s  = 


3  sinh  knL 

_ P  A _ 

2k  cosh2  kph  sinh3  k  h 
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3  p 
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and  cop  is  the  spectral  peak  angular  frequency  and  T  is  the  gamma  function.  It  should  be  noted 
that  the  analytical  results  of  Myrhaug  et  al.  (2009)  have  not  been  compared  to  experimental 
datasets  and  should  be  used  with  caution. 

VEGETATION  MORPHOLOGY:  As  discussed,  drag  coefficients  are  often  used  as  lumped 
calibration  parameters  accounting  for  different  vegetation  types,  material  properties,  and 
morphologies.  In  actuality,  individual  plants  respond  to  varying  environmental  conditions  and 
adjust  their  morphology  and  structure  accordingly.  Intra-  and  inter-species  competition  then 
governs  community  structure  and  composition.  These  individual,  population,  and  community 
level  processes  vary  over  space  and  time  to  result  in  a  patchy  landscape  of  variable  hydraulic 
roughness.  This  section  reviews  ongoing  mechanisms  and  processes  briefly  at  each  of  these 
scales  of  analysis  and  the  resulting  effects  on  wave  attenuation.  Complex  feedback  mechanisms 
also  exist  between  the  wave  environment  altering  vegetation;  however,  these  are  beyond  the 
scope  of  this  review  (See  Koch  et  al.  2006  for  review). 

At  the  individual  level,  vegetation  effects  on  hydrodynamics  are  the  result  of  plant  material 
properties,  surface  textures,  and  architecture.  The  primary  material  property  of  interest  is  the 
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rigidity  or  flexibility  of  the  plant  material,  often  reported  as  Young’s  modulus  of  elasticity 
(Ghisalberti  and  Nepf  2002,  Freeman  et  al.  2000).  Of  particular  interest  in  wave  environments  is 
the  material  density  or  specific  gravity,  which  determines  an  individual’s  buoyancy  and  location 
in  the  water  column  (Koch  et  al.  2006).  Skin  friction  of  plant  material  varies  with  leaf  and  stem 
conditions,  and  in  some  cases  exposed  roots  also  contribute  (e.g.,  mangrove  swamps).  Lastly, 
individuals  adjust  their  above-ground  and  below-ground  architecture  depending  on  ambient 
environmental  conditions  such  as  light  or  soil  density.  Above-ground  architecture  manifests 
through  changes  in  canopy  shape,  frontal  area,  plant  height,  and  many  other  measures  of  plant 
structure.  Although  rooting  conditions  do  not  influence  wave  attenuation  directly,  they  control  an 
individual’s  susceptibility  to  bending  and  washout  (Peltola  2006),  which  both  play  key  roles  in 
detennining  frictional  forces  (Fonseca  et  al.  1982). 

Plant  community  composition  and  structure  varies  widely  owing  to  spatially  and  temporally 
variable  environments,  competition  between  individuals  of  a  given  species,  competition  between 
species,  and  interactions  with  the  faunal  community  (e.g.,  herbivory).  Coastal  vegetation 
assemblages  vary  from  stands  of  near  monocultural  composition  (e.g.,  mangroves,  invasive 
phragmites)  to  diverse  communities  with  many  taxa  (i.e.,  groups  of  plants,  Figure  6).  Although 
community  composition  is  well-acknowledged  as  an  important  factor  influencing  wave 
attenuation,  to  date,  there  have  been  few  studies  examining  multiple  taxa  simultaneously. 
Moreover,  the  density  and  spatial  arrangement  of  the  community  has  been  demonstrated  to  be  an 
influential  driver  of  wave  attenuation  (Tanino  and  Nepf  2008,  Nepf  1999). 


Figure  6.  Two  marsh  plant  communities  in  Currituck  Sound,  NC:  (a)  a  diverse  assemblage  composed  of 
Juncus  romerianus,  Spartina  spp.,  and  numerous  other  taxa  and  (b)  a  near  monospecific 
stand  of  Phragmites  australis. 


In  addition  to  complex  individual,  population,  and  community  level  processes,  spatial  and 
temporal  variability  in  abiotic  (e.g.,  chemical  or  physical)  conditions  can  influence  vegetation 
morphology,  and  as  a  result,  its  effects  on  wave  attenuation.  In  a  flume  study  of  riverine  shrubs, 
Jarvela  (2002)  observed  that  leaf-on  conditions  could  produce  double  the  friction  force  of  leaf- 
off  conditions  for  the  same  species.  The  magnitude,  frequency,  duration,  and  timing  of 
disturbances  (e.g.,  erosion  events,  herbicide  application,  low  light  conditions)  can  also  lead  to 
changes  in  morphology  or  community  structure  such  as  vegetation  washout  or  scouring,  age 
structuring,  or  differential  growth  rates.  Table  3  provides  a  few  examples  of  individual, 
community,  and  spatio-temporal  processes  affecting  wave  attenuation  in  four  common  coastal 
ecosystem  types. 
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Table  3.  Examples  of  the  influence  of  individual,  community,  and  spatio-temporal 
dynamics  on  wave  attenuation  by  common  coastal  plant  communities. 

Coastal  Plant 
,  Community 

Individual  Properties 

Community  Properties 

Spatio-Temporal  Dynamics 

Macroalgae 
(e.g.,  kelp) 

Seaweed  fronds  becoming  more 
slender  as  they  increase  in  size 
(Gaylord  and  Denny  1997) 

Seaweed  fronds  “recline” 
against  each  other  to 
withstand  high  velocity 
(Denny  and  Gaylord  2002) 

Kelp  buoyancy  inducing  stem 
tension  at  high  depths  (Denny 
and  Gaylord  2002) 

Submerged  aquatic 

vegetation 

(e.g.,  eelgrass) 

Seagrass  bending  in  response  to 
increased  velocity  (Fonseca  et  al. 
1982) 

Natural  row  alignment  of 
seagrasses  in  response  to 
flow  directionality  (Fonseca 
et  al.  2007) 

Distribution  of  seagrass  beds  is 
often  due  to  hydrodynamic 
disturbance  (Koch  et  al.  2006) 

Wetland  or  Marsh 
(e.g.,  cordgrass) 

Differences  in  plant  morphology  of 
Spartina  alterniflora,  Juncus 

roemerianus,  and  Ditichlis  spicata 
result  in  different  velocity  profiles 
(Leonard  and  Luther  1995) 

Increased  bulrush  density 
reduced  wave  transmission 
in  both  field  and  laboratory 
tests  (Tschirky  2000) 

Wetland  vegetation  plays 

different  roles  in  wave  setup 
depending  upon  water  depth 
(Irish  et  al.  2008) 

Forested  swamp 
(e.g.,  mangrove) 

Exposed  rooting  zone  of 

mangroves  inducing  high  frictional 
resistance  at  low  water  levels 
(Massel  et  al.  1999) 

Mangrove  density  impacts 
the  role  of  wave  breaking  on 
overall  attenuation  (Vo- 
Luong  and  Massel  2008) 

Mangrove  swamps  create 

asymmetric  tidal  currents  and 
sustain  tidal  channels  (Wolanski 
1992) 

This  combination  of  individual,  community,  and  spatio-temporal  processes  alters  wave 
attenuation  through  two  primary  parameters:  area  exposed  to  flow  and  drag.  As  reviewed  by 
Fischenich  and  Dudley  (2000),  area  may  be  measured  through  a  variety  of  means  depending  on 
the  analytical  roughness  fonnulation  required.  As  shown  in  previous  sections,  drag  is  often 
lumped  together  analytically  because  of  the  numerous  confounding  variables  and  processes.  This 
CHETN  has  reviewed  drag  fonnulations  extensively  with  respect  to  wave  attenuation;  however, 
the  reader  is  encouraged  to  examine  the  numerous  formulations  from  riverine  literature  (e.g., 
Fischenich  2000,  Freeman  et  al.  2000,  Kees  et  al.  2010,  McKay  and  Fischenich  2011). 
Furthennore,  vegetative  roughness  in  rivers  has  been  accounted  for  historically  with  calibrated 
photographic  approaches  (e.g.,  pictorial  comparison  in  Bames  1967  and  Arcement  and  Schneider 
1989),  an  approach  not  yet  applied  to  coastal  environments. 

CONCLUSIONS:  The  ability  of  coastal  plants  to  dissipate  wave  energy  and  wave  heights  in 
low-energy  environments  is  demonstrated  and  documented  in  both  field  and  laboratory  studies. 
Wave  attenuation  is  highly  dependent  on  the  plant  biomechanical  and  wave  characteristics,  as 
evidenced  by  the  high  variability  in  wave  attenuation,  and  generalizing  vegetation-wave  behavior 
is  extremely  difficult.  Capturing  the  wave  transfonnations  induced  by  vegetation  is  critical  for 
predicting  coastal  hydrodynamics  accurately,  and  modeling  these  transfonnations  is  required  to 
quantify  the  benefits  of  vegetation  for  stonn  damage  reduction.  This  has  lead  to  the  development 
of  a  number  of  wave-vegetation  hydrodynamic  fonnulations,  which  account  for  the  influence  of 
plant  structure  on  wave  attenuation  through  a  calibration  variable,  often  CD,  that  is  adjusted  to 
minimize  the  difference  between  predicted  and  observed  wave  heights.  These  models,  and  the 
datasets  they  are  calibrated  with,  help  further  our  understanding  of  wave  attenuation  by  wetlands, 
and  serve  as  the  foremost  methods  in  describing  vegetative  dissipation  in  coastal  modeling  and 
planning.  However,  a  number  of  knowledge  gaps  still  exist: 


•  Quantifying  the  impact  of  stem  interactions  (based  on  spacing  and  structure/randomness  of 
plants) 

•  Validating  existing  dissipation  fonnulations  for  storm  conditions  (large  waves  and  surge) 

•  Quantifying  the  role  of  wave  period  on  dissipation  over  a  wide  range  of  depths  and  periods 
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•  Determining  drag  coefficients  based  on  wave  and  plant  parameters  instead  of  through 
calibration 

•  Including  seasonal  impacts  on  dissipation  potential 

•  Bridging  the  gap  between  the  larger  scale  at  which  wave  models  are  applied  (meters  to 
hundreds  of  meters)  and  the  much  smaller  scale  of  the  plant  features  and  dissipation 

•  Detennining  when  vegetation  lays  over,  breaks  off,  or  uproots  and  thus  becomes  ineffective 
to  dissipate  waves 

•  Quantifying  the  role  of  community  composition  (i.e.,  composition  and  configuration  of 
multiple  species  of  plants)  in  vegetative  roughness 

•  Validating  fonnulations  in  combined  wave  and  current  environments 

•  Detennining  the  interaction  of  wave  nonlinearity  and  dissipation  for  shallow  depths 

Existing  studies  of  wave  dissipation  by  vegetation  provide  a  range  of  analytical,  empirical,  and 
numerical  models,  but  current  methods  require  calibration  and  application  within  the  narrow 
range  of  available  lab  and  field  data.  Future  studies  need  to  expand  the  range  of  the  data  as  well 
as  generalize  model  fonnulations. 

SYMBOLS: 


A 

As 


Q-rms 

cD 

D 

Cfirms 

Cm 

Cp 

d 

Fd 

Ft 

F 

1  max 

Ft 

Fx 

F\ 

Fz 

h 

H 

H0 

Hrms 

k 

K 


k 

kA 

kp 

L 

Is 


cross-sectional  area  of  the  body  perpendicular  to  the  flow 

area  of  stem  normal  to  relative  velocity 

root-mean-square  wave  amplitude 

empirical  bulk  drag  coefficient 

drag  coefficient  for  a  singular  element 

root-mean-square  drag  coefficient 

inertial  coefficient 

empirical  plant  coefficient 

average  plant  stem  diameter 

drag  force 

inertial  force 

maximum  horizontal  drag  force  within  a  wave  cycle 
hydrodynamic  forcing 

horizontal  force  per  unit  volume  on  a  stem  array 
force  exerted  by  a  single  stem  on  the  fluid 

vertical  force  per  unit  volume  on  a  stem  array 

water  depth 

local  wave  height 

incident  wave  height 

root-mean-square  wave  height 

wave  number 

Keulegan-Carpenter  number 
exponential  damping  coefficient 
nonnalized  exponential  damping  coefficient 
wave  number  associated  with  spectral  peak 
wave  length 

average  plant  stem  length 
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N 

n 

P 

Q 

Red 

t 

T 

T 

1p 

tv 

U 

uc 

M  max 

Ur 

Uv 

V 

W 

X 

z 

r 

As 

c 

v 

P 

Pv 

T 

CO 

COp 

z 


plant  stem  density 

ratio  between  group  speed  and  celerity 
dynamic  pressure 

modified  Keulegan-Carpenter  number 

stem  Reynolds  number 

time 

wave  period 
peak  wave  period 
stem  width 

horizontal  wave  particle  velocity 

characteristic  velocity  acting  on  the  vegetation 

maximum  horizontal  velocity  within  a  wave  cycle 

relative  velocity  between  uv  and  u 

stem  swaying  velocity 

body’s  volume  perpendicular  to  the  flow 

vertical  wave  particle  velocity 

cross-shore  propagation  distance  through  a  vegetation  field 
vertical  coordinate 
gamma  function 
average  stem  spacing 

horizontal  displacement  of  the  vegetation  stand 
kinematic  viscosity 
fluid  density 

density  of  vegetation  element 
ensemble  interaction  coefficient 
wave  angular  frequency 
peak  spectral  angular  frequency 
vertical  coordinate 
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